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Abstract 

The Kraichnan rapid advection model is recast as the stochastic dynamics of tracer 
trajectories. This framework replaces the random fields with a small set of stochastic 
ordinary differential equations. Multiscaling of correlation functions arises naturally 
as a consequence of the geometry described by the evolution of N trajectories. Scaling 
exponents and scaling structures are interpreted as excited states of the evolution 
operator. The trajectories become nearly deterministic in high dimensions allowing for 
perturbation theory in this limit. We calculate perturbatively the anomalous exponent 
of the third and fourth order correlation functions. The fourth order result agrees with 
previous calculations. 

1 Introduction 

Although most researchers in the field of turbulence agree that the structure functions 
of the velocity field exhibit multiscaling, there is no well understood mechanism for 
this phenomenon which would allow, even in principle, a systematic calculation of the 
values of the multiscaling exponents Q. The Kraichnan rapid advection model || is 
a simplified model for turbulent advection of a passive scalar in which a mechanism 
for multiscaling has been identified. Although the relevance of this mechanism to 
multiscaling in fully developed turbulence in still unclear, the study of multiscaling in 
Kraichnan's model presents a new approach and therefore has attracted a lot of interest 
recently § g |, |. 

One of the main advances in studying this model was the identification of anoma- 
lous scaling contributions to correlation functions that are zero modes of certain linear 
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partial differential operators ||. Although this identification has allowed for the de- 
velopment of perturbation theories for the computation of the anomalous exponents, 
the meaning of these structures is, in our view, still somewhat obscure. The purpose 
of this paper is to reinterpret the anomalous exponents as eigenvalues of an evolution 
operator for the relative shape of trajectories of particles advected by the flow. The 
zero-modes arise naturally as the adjoints of the set of eigenfunctions of this evolution. 

Our method is to consider the evolution of N particles with the flow. These N 
particles define a configuration that has an overall scale and a normalized shape. Under 
the flow dynamics, the overall scale tends to increase while the shape will asymptotically 
be described by a stationary measure. We consider the rate at which a distribution 
of initial shapes tends towards its stationary distribution. It is important to consider 
this rate in terms of the overall scale. Doing this we find the decay rate towards the 
stationary state to be a combination of power laws in this scale. We identify these 
power laws with the anomalous exponents of the passive scalar. Thus we view the 
anomalous scaling as a dynamic property of the evolution of trajectories 0. 

This point of view treats the anomalous scaling as a geometric property of the 
trajectories of fluid particles in the flow, without explicit reference to the forcing. In 
fact, one does not need to define a passive scalar field since all the properties we are 
interested in are included in the shape dynamics. Note that this explains why the two 
point correlation functions do not have anomalous scaling — the geometry of a two 
point configuration is trivial, since it is completely described by the distance between 
the points. 

We will show that this picture formally contains all the ingredients of anomalous 
scaling in the dynamics of passive scalars. Furthermore, it is useful as a calculational 
scheme. We will demonstrate this by presenting an alternative derivation for perturba- 
tion theory in large dimensions. We believe that the physics of passive scalar advection 
in large dimensions becomes more transparent in this framework. We calculate the pre- 
viously unknown perturbative correction to the anomalous exponent of the third order 
correlation function, and recover the known result for the fourth order case ||. 

In section 2, we present the path integral formalism of passive scalar advection, and 
show how correlation functions can be written as averages over Lagrangian trajectories. 
In section 3 we exploit this picture to identify a mechanism for generation of anoma- 
lous scaling, and show how the anomalous exponents are related to eigenvalues of an 
evolution operator. In section 4 we apply this picture to calculate perturbatively the 
first correction to the anomalous exponents in large dimensions. Section 5 is devoted 
to summary and conclusions. 

2 Lagrangian path formulation of passive scalar advection 

The dynamics of a passive scalar 9 advected by an incompressible velocity field u are 
described by 

d t 9(r, t) = kV 2 9(t, t)-u- V0(r, t) + /(r, t) (1) 

where k is the molecular diffusivity, and / models the injection and extraction of scalar 
by external sources. The properties of u and / are given, presumably as the result of 
some physical process. 

In the rapid advection model j|] the velocity field is taken to be a random Gaussian 
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field, delta correlated in time and self similar correlations in space: 

<(u(r, t) - u(r',t)) (u(r, t') - u(r', t'))) = h(r - r')5(t - if), (2) 
where the "eddy-diffusivity" tensor h is defined by 

where the coefficients are chosen such that V • h = 0. The forcing / is also taken white 
in time, Gaussian, and is characterized by a single large scale L, 

{f(r,t)f(r',t')) = E(r-r')5(t-t'), (4) 

where the function H(r) is nearly constant for r <C L and decays rapidly for r > L. 

The rapid advection model is theoretically attractive due to the fact that the equa- 
tions for the correlation functions are linear partial differential equations as opposed to 
the hierarchy of nonlinear integro-differential equations for fully developed turbulence. 
Some theoretical progress has been made in the case of the rapid advection model and 
it is predicted that the 2nth order correlation function behaves as 

F 2n (ri, • • • , r 2n ) = L"( 2 -«)(c + • • • + c k (r/L)^F 2n (v) + •••), (5) 

where r <C L is a typical distance between points, and r denotes a set of dimensionless 
variables describing the configuration of the 2n points. According to the existing the- 
ory, the terms in the expansion (||) are of two types: Either they are inhomogeneous 
terms, whose scaling is determined by dimensional analysis, or they are the homo- 
geneous solutions of moment equations, whose scaling exponents are not constrained 
by dimensional arguments. The exponents and scaling functions are expected to be 
universal, but not the c coefficients, which depend on the details of forcing. 
The exponent ^ 2n which dominates when considering structure functions 

S 2n (ri,r 2 ) =< (0(n,t) - 9(r 2 ,t)) 2n > (6) 

is the scaling of the leading term in (|5|) which depends non-trivially on all of the In 
variables. Anomalous scaling occurs when this term is due to homogeneous solutions. 

As the path integral formalism for equation ([l]) is well documented H, we will not 
detail it here. The main conclusion is that one may write its solution formally 

9(r,t)= f dtf <f(r'(t'),t')> v (7) 
with the trajectory r' obeying 

r'(i) = r (8) 
d t ,r'{t') = 11(1/(0,0 +V^(0 (9) 

and r] is a vector of zero- mean independent Gaussian white random variables. 

The physical meaning of this equation in simply understood if we consider the k = 
case: the trajectories y are simply the characteristics of the partial differential equation 
(PDE) (H), namely the La grangian trajectories of the flow field u. If we want to know 
how much scalar is at point (r, t) we need to see what the value of 9 was where this 
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point came from, and add the accumulating forcing along the path traced backward 
in time Q. When k is positive the trajectories are randomly perturbed Lagrangian 
trajectories and the scalar 6 is obtained as an average over rj. This procedure generates 
the diffusion term in equation (|l|). 

Using the path integral formalism, we want to compute correlation function of the 
type 

F n (r 1 ,...,r n ) = (9(r 1 )---9(r n )) uJ , (10) 

where {A) u j denotes an average over realizations of velocity and forcing. Substituting 
each factor of #(rj) by its representation (Q) we can perform the / averages explicitly, 
yielding 

F 2n (ri, r 2n ,t) = (Jl^ dh ■ ■ ■ dtnlE^ih) - r'^O) • • • 



xH(r' 2n _ 1 (t n ) - r' 2n (i n )) + permutations 



while the odd moments vanish. Each of the trajectories obeys an equation of the 
form (^), where u as well as rj is a stochastic variable, whose correlation function is 
given by (||). The procedure just described is a Langevin dynamics of ./V particles, with 
correlated noise. Incompressibility of u implies that eq. (||) should be interpreted in 
the Ito sense. The corresponding Fokker-Planck equation with respect to both r] and u 
is precisely the well-known homogeneous Kraichnan PDE for the 2nth moment flOfl . 

The representation of F 2n given in eq ( |i~l"l) can be turned into a recursive expression 
by ordering the t\,---,t n integrations such that t\ > t 2 > ■ ■ ■ > t n . Since eq ( |Tl| ) 
contains all permutations in r' fc , this ordering can be done without changing the results. 
The times t 2 , ■ ■ ■ , t n can be integrated out to give 

F 2n (ri,..., r 2n ,t) = 

t r i 

(/-oo dt i[ H ( r i(*i) - r 2( t i)) i? 2n-2(r / 3(ti),---,r 2n (ti)) + permutations] ) U)V . 

This form for F 2n is useful because it only contains a single time integration, so that 
one needs to follow a joint trajectory of 2n particles at the same time and integrate 
along this trajectory the forcing defined by the terms in the square brackets in eq (|l2|). 

We intend to study the small k limit of the problem. If a given trajectory realiza- 
tion does not contain a 'near hit' among the trajectories, we expect that the k = 
trajectories are in some sense near the small k trajectories. The only case where k = 

differs radically from the small k limit is when two (or more) trajectories become closer 
i 

than 0(k 2 -z). As we are interested only in separated correlation functions, we neglect 
the corrections to the k = limit in what follows. 

All the existing studies of the Kraichnan model are based on analyzing the Kraich- 
nan linear PDEs. Most of these studies identify the anomalous contributions as scale 
invariant functions of the coordinates which are annihilated by the Kraichnan operator, 
thus they were termed zero modes. Most researchers accept the picture that correlation 
functions are obtained as a sum of zero modes with different anomalous scaling expo- 
nents and a contribution which is directly related to the forcing which is not anomalous. 
A notable exception to this approach is the Kraichnan closure, based on an appealing 
but unproven conjecture on conditional expectations of diffusive moments Q. 

The common feature for these approaches is that the passive scalar is treated as 
a field. Here, we propose an alternative viewpoint based on eq. (Ill]), in which all the 
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properties of low order moments are extracted by looking at a small number of corre- 
lated stochastic trajectories. The representation (pd|) expresses F2 n as the expectation 
value of the forcing correlation function accumulated over the Lagrangian trajectories. 
Since the forcing drops off sharply beyond L these objects measure the correlation of 
the times for which the trajectory positions are composed of pairs of points which are 
within distance L. In particular the two point moment F2 is the average time during 
which two trajectories which end at given distance of each other stay within distance 
L. Excluding the case d = 2, £ = the trajectories will leave the vicinity of each other 
with probability 1, and the integrals in (|ll]) converge for each realization. 

It should be noted that the representation described above can be used in principle 
as a numerical scheme for computing directly the correlation functions F<m using a 



Monte-Carlo method without the need to generate the whole velocity field [11 



In what follows we show how one can eliminate the forcing and the outer scale 
from the problem. Anomalous scaling arises from a relaxation process, in which the 
configuration tends towards an asymptotic stationary distribution, while the overall 
scale increases. The outer scale serves as a reference point for this process, but is not 
intrinsic to the process itself. 



3 Shape evolution operator and its relation to multiscal- 
ing 

A component of the evolution of an initial configuration is a rescaling of all the co- 
ordinates which all increase on the average like t 1 ^ 2 ; this rescaling is analogous to 
Richardson diffusion for the present case The dynamic exponent £2 = 2 — £ is also 
the characteristic exponent of the second order structure function ||] . After factoring 
this overall expansion we are left with a normalized 'shape'. It is the evolution of this 
shape that determines the anomalous scaling. 

Consider an initial shape Zo with an overall scale so, that is, the configuration 
coordinates are Xo = sqLq. Note that the shape coordinates are constrained. The 
overall scale is a homogeneous function of degree one of the coordinates so = 5(Xo) 
but its exact functional form is to a large degree arbitrary. 

The shape evolves in time. We fix a scale s > so and examine the shape when 
the configuration reaches the scale s for the first time (we know that this occurs in 
finite time with probability one). Since the trajectories are random the new shape Z 
is taken from a distribution 7(Zo, Z, ^). We used scale invariance to deduce that the 
distribution is a function of — . 

For a shape picked from an initial distribution of shapes p so (Zo) with scale sq, the 
final shape distribution at scale s is 

p a (Z) = [ dZ 0/ 9 (Zo)7(Zo,Z,-) (13) 

For very large values of ^ we expect the final shape distribution to approach an 
asymptotic distribution /3q(Z). This distribution is invariant under the transformation 
(|l~3| ) and is an eigenfunction of 7 with eigenvalue 1. In general we expect a spectrum 
of eigenfunctions f3 n for 7 with eigenvalues a n (^). By successive applications of 7 it 
follows that a n (— ) = (— )~ A ™. Furthermore we assume that the eigenfunctions /3 n (Z) 
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form a complete set. Thus, any initial distribution of shapes p S0 (Z) can be written as 

Pso (Z)=5>„/3„(Z) (14) 



this distribution evolves into 



s x A " 



ft(Z) = V4 - Pm (z) . (15) 

Vso/ 

Consider an average of a correlation function F n (X) over p SQ , 

(F n \s ) = J dZF n (sZ)p S0 (Z). (16) 



We now use the representation (|12|) for F n (so, Zo). We split the integration along the 
Lagrangian trajectory at the first point where the overall scale reaches the value s. It 
follows from the shape evolution equation (]i~3| ) that we can write 

F n (s , Z ) = /(s , s, Z ) + y dZ 7 (z , Z, ^ F n (s, Z) . (17) 

/ is the average of the contribution a specific forcing accumulated along the path from 
so to s. The forcing is that of eq (|l2] ) so that for s small enough its L dependence is 
only through F n _2- Averaging eq. ( |i~7| ) with / o SQ (Zo) gives 

/ dZ[p so (Z)F n (s Z)-p s (Z)F n (sZ)) = (l( So , s ))=J2A m (Ft\s )-(-) Xm Ft ) (s)) 
J m \s J 

(18) 

Where 

= J p m (Z)F n (sZ)dZ (19) 

and 

(I(s ,s)) = J Pso (Z )I(s ,s,Z )dZ . (20) 

The average {I(so,s)} is dependent on L only through the behavior -F n _2 with L. 
Thus, this term cannot have anomalous 2'( 2_ ^)~^ n ) scaling. Since the A^s are arbi- 
trary, each term in the sum of equation (18) must be without anomalous 2 ( 2_ 0~Cn) 



scaling as well. The only way for this to occur is if the anomalous L-dependent part of 
Fn m \s)) is proportional to f(L)s Xm . Dimensional considerations imply that f(L) also 
has a power-law behavior that fixes the dimension of F n . 

Thus, the projection of F n over (3 m selects a particular anomalous component. 
We conclude that the exponents A m , defined as the decay rates toward the invariant 
measure, are precisely the anomalous exponents of F n . 

We can use the eigenfunctions (3 n to expand the operator 7 



s { s \ 171 

7 (Z ,Z,-)=E - /UZWZq). (21) 



7 is a non-Hermitian operator, and therefore we can expect only that (3 m and p m are 
biorthogonal families of functions. This means that if we expand F n in terms of the 
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functions /U m , a projection on (5 m will extract a single term of the sum, which has a 
pure anomalous scaling exponent. This leads us to identify the functions [x m with the 
zero modes of the Kraichnan operator, which were previously identified in || ||| as the 
anomalous scaling structures. 

We will use this formalism in the next section in an explicit perturbative calculation 
in large dimensions. This calculation demonstrates all the concepts discussed in this 
section, and shows how they can be used as a calculational method. 

4 Application: large d perturbation theory 

The simplicity of the large dimension limit is due to the following reason: random 
trajectories tend to separate faster the larger the dimension. This is because there are 
more transversal directions in which to diffuse. This property which is well known for 
Brownian motion, is also true for our case of correlated, Markovian paths. Therefore, 
in the limit of infinite dimensions the motion becomes a deterministic growth of the 
distance between each pair of points independently (note that the infinite number of di- 
rections available makes such a motion possible) . This motion is simply the Richardson 
diffusion phenomenon, referred to above, in which relative distances I increase like 

l~ (z + t )V(2-O. (22) 

This behaviour means that any initial iV-point configuration will evolve towards a 
regular N — 1 dimensional simplex, where all the relative distances are equal. 

When d is large but finite, trajectories fluctuate around the deterministic infinite 
dimensional limit (]22]), but the fluctuations are small, of order 0(l/y/d). This forms 
the basis of perturbation theory in 1 jd. 

Our procedure will be to examine the expectation value of a non-dimensional sym- 
metric combination of separations, which we denote a. The relaxation of (er) towards 
its asymptotic value is described by the set of relaxation exponents A m , see (^) 

(a\s) = J dZa(Z) 7 (z , Z, £) = £ (£) ™ a mf x m (Z ), (23) 

where Zq is the initial shape, implying an initial value for <r(Zo). The cr m 's defined as 

J a m = dZa(Z)/3 m {Z), (24) 

are numbers which only depend on the precise definition of a. We see from (^) that 
the functional dependence of {cr\s) on the initial shape Zq is given by the zero modes 
fj, m . An asymptotic expansion in large will thus provide us with the leading zero 
modes, and their scaling exponents. 

A term in the asymptotic expansion (cr\s) is expected to be, in the large d limit, of 
the form 

[ M (°)(Zo) + -^ 1 )(Zo) + •••]. (25) 

Expanding the exponent gives 

/ <■ \ ~ A(0) 1 <> 

[ M (°)(Z ) + -(-A« log(-) + /^(Zo)) + •••]. (26) 
\sqJ d s 
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Hence, the correction to the scaling exponent can be read off by looking at the loga- 
rithmic term in ^. Although nominally our expansion is in -7= it turns out that the 
first term vanishes (see below). 

In principle, the statement that a logarithmic contribution is the first term in an 
expansion of a power needs to be justified by showing that all higher order terms of 
order ^log(-^) n are consistent with the first term. In our case this is not necessary. 
Rather, we consider a perturbative expansion for the equivalent Fokker-Planck equa- 
tion. In the Fokker-Planck equation the eigenvalues play the role of the exponents X m . 
Since perturbation theory finds corrections in the eigenvalues, we conclude that the 
logarithms must sum to a pure power as expected. 

The perturbative analysis is going to carried out in terms of time as an independent 
variable rather than s. We thus need at some stage to transform to the s variable in 
which our theory is formulated. Since the dynamics are nearly deterministic, for a 
given s, t is very narrowly distributed around its Richardson diffusion value t ~ s 2 ~^. 
We show that to the order that was retained in the perturbation calculations it is 
sufficient to replace t by its Richardson diffusion value. 



4.1 Representation in terms of separations, and perturbation theory 

The basic variables we use for analysis in high dimensions are the n s = N(N — l)/2 
inter-point square separations q nm = (r n — r m ) 2 . Their time evolution is obtained from 
(|7|) by the rules of the Ito calculus 

d(jnm{t) 1/ \l\ , <->/ \ / \ n?nmW /j £ \ . ~ r,\ 

K-u m ) ) + 2(r n -r m )-(u n -u m ) = 2 (d-- — — — ) + % m (£), 



,. \ V ~-7fc —/ftp I / 1 — V It, ill! \ —ft — fft / f)C \ j -1 , >- , 

at \ ' t> a — 1 + 4 

(27) 

where the second equality serves to define the noise term rj nm , a zero mean, 5-correlated 
in time Gaussian process. The equations are again to be interpreted as Ito SDEs. The 
n s separations are subject to triangle inequalities which are preserved by the dynamics 

& 

The deterministic part in eqs. ([27]) gives simple Richardson diffusion with the expo- 
nent 2/(2 — £) = 2/^2, whereas all the non-trivial behavior is contained in the fj terms 
which couple between the different separations. When d is large the noise term in ( |27| ) 
becomes small with respect to the deterministic term, and this phenomenon serves as 
the starting point for perturbation theory. In order to demonstrate this it is convenient 
to factor out the super-diffusive behavior from q nm , defining 

Hum (*) = -ir-d{t + (t), (28) 



where r nm is a constant determined by initial conditions. Substituting in eq. (p?|) yields 
(here and in the following we neglect the higher order term — 1+0 coming from 
incompressibility) 



~ (4!m ~ Qnm) + -^=r] nm (t), (29) 



mm „ c , \inm mm i 

I — k, t + T nm v / y 1 d 

where 

\ -2/(2-0 

Vnm = ( — d(t + T nm )J rj nm . (30) 
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The covariance of Tj nm is (using the abbreviation q nm = (t + r nm ) 2 ^ 2 ^>q nm ) 
(Vnm(t)Vkl(t')) = ir(t + r nm )- 2/ ^{t + r n ,)~ 2/C2 

x (q n l ~ Qml ~ Qnk + Qmk)^ 2 ~ qJ^ ~ + ~ • 

(31) 

This covariance is of order 1 in d, so that the stochastic term in eq. ( |29| ) is indeed small 
for d S> 1. Thus, the dynamics of eq. (29) are dominated by the attractive fixed point 
Qnm = 1 of the deterministic term, and the random term causes small fluctuations 
around it. 

Perturbation theory is performed in a straightforward manner by expanding 

Qnm = 1 + —F=^Qnm + ~^Qrirn + • • • 1 (32) 

and ^ 

^ nm = Vd^" 1 ~*~ d ^ ^ nm ^ kl + • • • • (33) 

The noise terms are white Gaussian zero-mean processes. The correlation matrix of 

rjnm, is obtained from ( f3l| ) by substituting s nm — ► 1 — <5 nm . The correlations involv- 

(2) 

ing r]nm are more complicated, but are not needed to the order that we keep in the 
calculations. 

The terms in the perturbation series (p2|) obey linear inhomogeneous equations: 



c,W = i M) , „(i) (-34) 

'inm 1 1 'dnm 1 Inm J v u ^v 
" ~r Trim 

and ^ 

<?nm = ~~ 7 , (^nm + T^im) 2 ) + X] ^nm,kl^kl ' (35) 
1 ~r "nm * ^ 



with respective solutions 



= — T ^ + T«m)^(0 , (36) 
t + T nm JO 



and 



= 7—^ / dt' 

* ~r 'nm J 



■|gW ( t )2 + {t' + r nm ) £ uS^OsSV) 



(37) 



4.2 Results for three-and four-point functions 



In this section we are going to use the results of perturbation theory derived above to 
demonstrate multiscaling in the three- and four-point correlation functions. For this 
purpose we make the following (arbitrary) choice 

a = n a |r wm - g "V (38) 

{l^nm Qnm) 

The definition of the overall scale is chosen to be 



Qnm ■ (39) 



N 

n<m 
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It is shown in the following subsection that replacing t by s& is correct to the order 
which we keep in the following calculation. 

Substituting the perturbation series (]3^) to order 1/d gives the expansion a = 
+ + + ' ' ■' w here we have defined 



r (0) 



n. 



l^nm Pnm 

(T o 2/C2 ) 



2 ' 



(40) 



a 



(l) 



a 



(0) 



Y- VC2 (1) v-> 2/C2 (1) 

2 2^ ram P nm Qnm _ 2^nm Pnm Qnm 



E 



4/C2 

nm Pnm 



E 



2/C2 

nm ^nm 



(41) 



(T 



(2) 



(o) 



l^nm Pnm (qnm) 



4/C2 
nm Pnm 

2/C2„(l) 



+3 



E 

2—inm Pnm Qnm 



E 



+ 2 



'v- 4/C2 (2) ^ 2/C2 (2) * 
Z^nm Pnm (/nm Z^ram Pnm Qnm 



Lnm Pnm 

4/C2 (1) 



Z^nm Pnm 

2/C2 fl) 



2-~i nm Pnm Qnm 2— mm Pnm Qnm 



2/C2 
nm P«m 



E 



4/C2 

nm Pnm 



E 



2/C2 

nm P«™ 



(42) 



(We used the abbreviation p nm = t + r nm .) 

When i 3> r nm , (cr) is dominated by contributions from the leading zero modes, 
and we therefore check its asymptotics in this limit. In the limiting case d = 00 we 
have, for t — > 00 



Enm 1 + (4/C2) (T„ m /t) + (2/C2)(4/C2 " 1)(W*) S 



(Enm 1 + (2/C2)(T„ m /t) + (1/C2)(2/C2 " 1)(W*) 2 ) 2 
-1 4 v -» / Tram. "N I 

nrr 



KC2) 5 



E 



t 



(43) 



In the special case N = 3, n s = 3 eq. (43) can be written as 



a 



(0) 



1 + 



_8_^ 
9C| t 2 



(44) 



where Z3 is the leading three-point zero mode, 

■^3 = 2-s ( T nm 

m < I 



Tnlf 



(45) 



is thus a linear combination of zero modes (up to order t 2 ) as expected. 

The 0(1/ y/d) contribution vanishes since (rf^j = 0, so that we proceed to 

examine the 0(1/ d) term. Up to order i~ 2 , (a) contains, in addition to the constant 
term, a term which is proportional to logt/t 2 . This term reflects the dependence of 
the scaling exponent on d, and we will therefore display only this term explicitly. 

The first step consists of calculating moments of the s functions: 



„(1) 

Hnm 



(t + Tn 



dt'dt"(t' + r nm )(t" + r nm ) (4l{t)^l{t' 



(46) 



C2 



(t + T n 
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T 2 
' nm 



2 

nm y 



(47) 

(The term proportional to rj^ drops in averaging due to the Ito convention). Thus, 
terms involving only one separation do not contribute logarithmic terms to this order. 

The large t asymptotics of the cross correlation of two separations with one common 
vertex is 

2 log* ( 2 A [ 1 (Jl t _ 2 > 1 , ( 2 A (ji , _2 _2 > ( 2 



> \ r nm + r nz) + ~~ 1J ( T nm + r ni ~~ r mz) + ~~ 1 ) T nmJnl 



C2 /4 \ 
Summing over permutations gives 



(48) 



^-'A5 T 'r' («) 

We conclude that the only terms in (&^/ which generate logarithms are the ones 

involving cross correlations of q variables, i.e., the last two in eq. (fi^). Of these terms 
we need only keep the leading order in 1/t, giving 

(2)\ , ci 4^2 ^ Z' 2 , ^ „ log* 



Combining the contributions from (44) and fl50| ) yields 

(a) ~ 1 + 9cf^ (1 -2( v 2^ + 1 J^r ) - (51) 



which implies, after substituting t —* s 2 ^ 



C 3 = 2(2-0 + !j(4-0 + O(-^) . (r.2) 



2d v sy Vd 3/2 

A similar analysis is needed for the case N = 4, n s = 6. The infinite dimensional 
limit becomes in this case 

■^-^W^- (53) 

where the leading order 4-point zero mode is 

Zi = ( T «™ ~~ T ki) 2 ■ ( 54 ) 

n<Cm, k-Cl 

We also need cross correlations of separation variables without a common vertex, whose 
large t asymptotics are 



(2\ 2 (2 Alogt. , 2 , Ktn 

VC2/ VC2 ~~ / ~ Tml ~ Tnk + T ' 
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and summing over permutations gives 

2 



£ (|) (§-l)^<-2% + *). (56) 



n<m, fe<Z 



As in the 3-point case we collect the contributions from (|53|), (^), and ( p6|) giving 
1 1 /_ _ log* 



w ~ 1 + 9ci?l Z3+z ' + i"« 

and the dependence on Z4 implies that 



4 + + 



(57) 



Ci = 2(2-0-f + 0(^), (58) 

in agreement with . Eq. ([5?]) yields further information by checking the contributions 
proportional to Z3; the relative strength of the logarithmic terms is different from 
that obtained in the three point analysis, see (^Tj). The reason is that generates 
contributions proportional to Z3 at first order, so that it is no longer an approximate 
zero mode when d is finite, and the degeneracy between the 3-point and 4-point zero 
modes breaks down. It is not difficult to build the 4-point zero mode Z4 which continues 
correctly from infinite to finite dimensions as a linear combination of Z3 and Z4 

Zi = Z 4 + -^-Z 3 . (59) 
4.3 The t —y s transformation 

The previous analysis relied on calculations as a function of time, and the transfor- 
mation to representation as a function of total scale was carried out by the simple 
substitution t — > s^ 2 . The purpose of this subsection is to show that corrections due to 
non-trivial dependence of s on t do not contribute terms of the order which was kept 
above, i.e. 0(1 /t 2 ) and 0(log t/(dt 2 )). 

It should be emphasized that even when d is very large, s is not a monotonic 
function of t, and this is evident even in first order perturbation theory. However, the 
size of the time interval in which one is likely to find the same of value of s becomes 
very small. In order to estimate the width of this interval, suppose s(t) = s. We ask 
for which values of r it is likely that s(t + r) = s as well. We know that for small r 

s(t + r) ~s + Ct + ^^lu, (60) 
yd 

where u is a random variable with O(l) variance, and B and C are some order 1 num- 
bers. The deterministic part is monotonic and may be compensated by the fluctuating 
part to give a solution only if r = 0(l/d). This is our estimate for the width of the 
multiple solution region. Replacing the earliest solution for t(s) by some other arbi- 
trary solution thus induces an 0(l/d) error in t, which is too small to affect the results 
to the order that we keep. 

In addition, up to 0(l/d) we may assume that s(t) is a one-to-one mapping, and 
may invert the relation to obtain t(s). s may be written (up to a constant multiplier) 
as 

s(t) ~ £(t + r rem )2/C 2 (l + -^qffi(t) + • • •) . (61) 
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Inverting to express t clS EL function of s gives 



t = t (s) [ 1 - ttTT" V T nm -= V q£ m (t (s)) ) + higher order terms, (62) 

V *0C s ; C2 n m Vd nm / 

where to = s& . We see that to 0(1/ Vd), t may indeed be considered as a single valued 
function of s which is however randomly shifted with respect to the zero order estimate. 
We are now going to show that this correction does not contribute to the perturbative 
results either. 

We substitute t(s) in the asymptotic expression for (see eq. 43, keeping the 
dependence on r nm implicit), 

- (0) ~ 1+%^ ~ 1 + f 1 + ~ E ^ m + E «a(*d) + ^(E e^0)) 2 ) 

J V H> C,2 „ m Va nm " nm / 

(63) 

It follows after averaging that corrections to to wu l contribute only terms of 0(l/tg) 
and O(logio/(*o))- 

(<7^ ^ , which vanishes when one takes t = to as in the previous subsection, is non- 
zero when corrections to to are taken into account. The leading non-zero terms are 
proportional to 

1 M 

1 \ ^ T nm q n m \ - (1) , \ 

-75 E — : E %i (*o) . ( 64 ) 

whose expectation is of order O(l/(s/dto)) and 0(<i~ 3//2 logtoAo)- Finally, it is clear 
that corrections to to m ( <T< ' 2 ^) cann °t change the the leading logarithmic behavior. 
We conclude that to the order to which the above calculations we carried out it is safe 
to take t = to = s& . 



5 Conclusions and further applications 

In this paper we have presented a path oriented approach to the study of passive 
scalar advection. The main difference between this approach and other studies is 
that we do not need to consider the whole flow (and scalar) field, rather, an iV-point 
correlation function is described via the evolution of N Lagrangian trajectories. This 
allows us to study ordinary differential equations (albeit stochastic ones) instead of 
partial differential operators. The anomalous scaling is due to relaxation towards 
an invariant distribution of the instantaneous shapes. Since the relaxation rates are 
associated with excited states of an evolution operator these scaling exponents are not 
related to the normal scaling and dimensional reasoning cannot be applied. 

The Lagrangian trajectories become nearly deterministic in very large dimensions 
(up to random solid body rotations). This property explains why the dynamics simpli- 
fies in this limit, and also serves as a starting point for an expansion. Using perturba- 
tion theory all the ingredients of the new anomalous scaling picture were demonstrated 
in an explicit and concrete manner, and perturbative corrections to the anomalous 
exponents were calculated. 

In addition to the perturbative application it is possible to use the same concepts for 
numerical Monte-Carlo simulations. The main numerical task is to generate the random 
trajectories, which can be done using standard methods. However, in preliminary 
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studies severe problems of convergence prevented us from obtaining precise results, 
and therefore the presentation of the numerical application is postponed |11, 12]. 



It would also be interesting to see how the other integrable limits of this model 
(£ — ► and £ — ► 2) appear using the present ideas. 

We would like to thank Uriel Frisch, Zeev Olami, Itamar Procaccia and Massimo 
Vergassola for discussions. O.G. thanks Krzystof Gawedzki and Itamar Procaccia, the 
organizers of the Turbulence Workshop in I.H.E.S in April 1997, where many of the 
ideas in this paper were developed. Most of this work was done while R.Z. was a CNRS 
post-doctoral fellow in the Laboratoire de Physique Statistique at the Ecole Normale 
Superieure (Paris). The Laboratoire de Physique Statistique is Laboratoire associe aux 
Universites Pierre et Marie Curie, Denis Diderot et au CNRS. 

References 

[1] U. Frisch, Turbulence, Cambridge University Press (1995). 

[2] R. H. Kraichnan, Phys. Fluids 11 945 (1968) 

[3] R. H. Kraichnan, Phys. Rev. Lett., 72 1016 (1994) 

[4] K. Gawedzki and A. Kupiainen, Phys. Rev. Lett., 75, 3834, (1995). 

[5] M.Chertkov, G. Falkovich, I. Kolokolov and V. Lebedev, Phys. Rev. E 52 4924, 
(1995). 

[6] O. Gat, V.S. L'vov, and I. Procaccia, Phys. Rev. E56, 406 (1997). 

[7] A similar point of view is presented in D. Bernard, K. Gawedzki, A. Kupiainen, 
e-print |omd-mat/970603"5| (1997). 



[8] See for example, I. T. Drummond, J. Fluid Mech. 123, 59 (1982), M. Chertkov, 
Phys Rev. E 55, 2722, (1997), P. M. Ginanneschi, e-print no. |chao-dyn/9704017 
(1997). 

[9] B. Derrida, R. Zeitak, Phys. Rev. E. 54 pp 2513 (1996). 
[10] R. H. Kraichnan, V. Yakhot and S. Chen, Phys. Rev. Lett., 75, 240 (1996). 
[11] O. Gat, I. Procaccia, and R. Zeitak, unpublished (1997). 

[12] The idea of using paths for passive scalar Monte-Carlo simulation was developed 
independently by U. Frisch and M. Vergassola. 



14 



